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Abstract. Models of black hole accretion disks are reviewed, with an 
emphasis on the theory of hard X-ray production. The following models 
are considered: i) standard, ii) super-critical, hi) two-temperature, and 
iv) disk+corona. New developments have recently been made in hydro- 
dynamical models of accretion and in phenomenological radiative models 
fitting the observed X-ray spectra. Attempts to unify the two approaches 
are only partly successful. 



1. Introduction 

Accretion disks around black holes are efficient radiators which can convert a 
sizable fraction of gas rest mass energy into radiation. The disk forms if the 
specific angular momentum of the accreting gas j ^> r g c, where r g = 2GM/c 2 
is the gravitational radius of the black hole. Then the radial infall is stopped 
by the centrifugal potential barrier at a radius ~ j 2 /r g c 2 , the gas cools and 
forms a ring rotating with Keplerian velocity. Further accretion is possible only 
if some mechanism redistributes angular momentum, allowing gas to spiral to- 
ward the black hole. A plausible mechanism is provided by MHD turbulence 
that develops due to the differential character of Keplerian rotation (see Balbus 
& Hawley 1998 for a review). Turbulent pulsations force each gas element to 
diffuse from one Keplerian orbit to another, and a ring of gas spreads out to 
form an extended disk. At the inner edge of the disk gas is absorbed by the 
black hole. The diffusion problem with the absorbing inner boundary has sta- 
tionary solutions, each specified by a net accretion rate, M. The redistribution 
of angular momentum is accompanied by a viscous heating. As a result, the 
binding energy of spiraling gas is dissipated into heat and can be radiated. 

In this review, we concentrate on steady disk models and their applications 
to active galactic nuclei (AGNs) and galactic black hole candidates (GBHs). 
Observed spectra of the black hole sources indicate that accretion disks usually 
do not radiate as a blackbody. In particular, a huge luminosity comes out in hard 
X-rays. The origin of the hard X-ray component has been in the focus of interest 
since the component was detected, and it remains as a challenge for theoretical 
models of accretion disks. We start with the standard disk model (Section 2) 
and then discuss the concept of advective disks (Section 3). In Section 4, we 
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discuss super-critical disk models. In Section 5, two-temperature models are 
reviewed. Finally, in Section 6, we summarize new developments concerning the 
disk-corona model. 



2. The standard model 

The standard model (see Pringle 1981 for a review) provides a commonly used 
simple description of an accretion disk. The model considers a thin gaseous 
disk of surface density S and density p = £/2iJ, where H is the disk height- 
scale. The disk rotates around the black hole with a Keplerian angular velocity, 
= (GM/r 3 ) 1 / 2 . H is regulated by the pressure in the disk, p, which supports 
the gas against the vertical component of gravity. The vertical balance reads 

r vk 

where c 2 = p/p and vk = &K r - The matter in the disk gradually drifts inward 
due to a viscous stress t rip = vpr{dSlYi/dr) where v is a kinematic viscosity 
coefficient. The stress is a fraction of the pressure, t rifi = ap, where a < 1 is 
assumed to be constant (Shakura 1972). Equivalently, v = (2/3)ac s H. The 
inward drift is governed by the angular momentum equation, 

— (n K r 2 ) = — (T rv r 2 ), T rip = 2Ht r(p = -vW K . (2) 

Here v is the accretion velocity. In a steady disk, 

27rruS = M = const, (3) 
which allows the angular momentum equation to be integrated, 

The constant C is determined by an inner boundary condition. In the very 
vicinity of the black hole, Keplerian rotation becomes unstable and gas falls into 
the black hole with a constant angular momentum. The radius of the marginally 
stable Keplerian orbit is r ms = 3r g for a Schwarzschild black hole. The standard 
model treats the transition to free fall by placing the inner disk boundary at 
fin = ?"ms and by assuming t r(p = at r = r Ui . Then C = ^/3r 9 . 
The accretion velocity may be expressed from (3) and (4) as 

» = ?V = «i S -, SE i_(iA 1/2 . ( 5 ) 

2r v K \ r J v ' 



The viscous heating rate (per unit area of the disk) is 

+ dft K 3 3M 2 

F = T rip r— — = -T rip {l K = — — U K b. (6) 



The standard model assumes that the bulk of the dissipated energy is radiated 
away locally, i.e., the local heating rate, F + , is balanced by the radiative losses 
from the two faces of the disk, 2F = F~ . The disk therefore remains cold and 
thin, c s <C vk, H <C r, and accretion is slow, v <C vk- 

As the binding energy is radiated away, the radiative efficiency of accretion 
is rj = bi n /c 2 , where 6j n is the specific binding energy at the inner edge. In the 
simplest version (eqs. [1-6]), a Newtonian gravitational field is assumed, with the 
potential GM/r. Then 6 in = GM/2r in = c 2 /12. A relativistic version (Novikov 
& Thorne 1973; Page & Thorne 1974) yields 0.06 < rj < 0.42 depending on 
the black hole spin. The highest efficiency is achieved for a maximally rotating 
black hole. Then the stable Keplerian disk spreads deep into the potential 
well, r ms — > r g /2, and this results in a high binding energy at the inner edge, 
b iQ = c 2 (l - 1/^3) « 0.42c 2 . 

The set of disk structure equations (1), (3), (5), and (6) should be completed 
by an equation of state and an equation specifying a mechanism of radiative 
cooling. In particular, if (i) the pressure in the disk is dominated by ionized gas 
with a temperature T, so that c 2 f» 2kT/m p , and (ii) the produced radiation is 
near thermodynamic equilibrium with the gas and gradually diffuses out of the 
disk, then F « ca r T 4 /V3r = {ca r m 4 /16V3k 4 )c 8 s /r, where a r = 7.56 x 10~ 15 is 
the radiation constant and 2r = <j<iYi/m p is the Thomson optical depth across 
the disk. It completes the set of equations and the disk parameters can be 
expressed as analytical functions of radius, each model being specified by M and 
a (Shakura & Sunyaev 1973). We will hereafter express M in a dimensionless 
form, m = Mc 2 /Le, where Le = 2irr g m p c 3 /a^ is the Eddington luminosity. 

When m exceeds ~ O.l(aM/M ) -1 / 8 there appears a region in the disk 
where the radiation pressure dominates the gas pressure. Then the disk thickness 
is determined by the vertical radiation flux, GMm p H/r 3 = Fa^/c, and the 
set of equations may be closed by the condition H = (3/4:)mr g S. The model 
then becomes viscously (Lightman & Eardley 1974) and thermally (Shakura & 
Sunyaev 1976) unstable. The instabilities are reviewed by Julian Krolik in these 
proceedings. 

Assuming that the disk radiates as a black body, F = aT 4 , one gets 

A maximum surface temperature is T* ~ 10 5 K for AGNs (M ~ 10 7 - 1O 9 M ) 
and T* ~ 10 7 K for GBHs (M ~ 1OM ). It is the major success of the model 
that T* approximately corresponds to observed UV emission in AGNs and X-ray 
emission in GBHs. The hard X-ray component would, however, correspond to 
a temperature ~ 10 9 K. The standard model may yield such high temperatures 
only if the accretion rate is near the critical value, rh cr = r? _1 , that corresponds 
to the critical Eddington luminosity. Near-critical disks can strongly deviate 
from the blackbody state provided a is large (Shakura & Sunyaev 1973). An 
overheating then occurs in the inner region of the disk where the inflow time- 
scale is shorter than the time-scale for relaxation to thermodynamic equilibrium. 
This effect is discussed in more detail in Section 4. 
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3. Advective disks: basic equations 

The key assumption of the standard model, that the dissipated energy is radiated 
away locally, is relaxed in the presently popular models of advection-dominated 
accretion flows (ADAFs, see Narayan, Mahadevan & Quataert 1998 for a re- 
view). There are two types of ADAFs: super-critical (see Section 4) and two- 
temperature (see Section 5.2). In both models, a large fraction of the released 
energy is stored in the gas and advected into the black hole instead of being 
radiated away. Advective disks are geometrically thick, and a two-dimensional 
(2D) approach would be more adequate to the problem. Most of the models are, 
however, based on the vertically integrated (ID) equations due to their simplic- 
ity. The a-parametrization of viscosity is used in the same manner as in the 
standard model, sometimes with modifications suppressing a at the sonic radius 
to avoid the causality problem (Narayan 1992). By solving the ID equations, 
one gets typical parameters of the accreting g clS clS cl function of radius. 

Advection implies a non-local character of the flow, so that one deals with 
a set of non-linear differential equations with boundary conditions. In the in- 
nermost region, the flow passes through a sonic radius, r s , where the radial 
velocity v exceeds the sound speed. A regularity condition must be fulfilled for 
any steady transonic flow at r = r s (cf. Landau & Lifshitz 1987). This condition 
reads N = D = where N and D are the numerator and denominator in the 
explicit expression for the derivative dv/dr. The regularity condition gives two 
extra equations that determine r s and impose a connection between M and the 
flow parameters at an external boundary. 

The transonic character of accretion is related to a special feature of a black 
hole gravitational field: the effective potential for radial motion with a given an- 
gular momentum has a maximum (see Misner, Thorne & Wheeler 1973). After 
passing the potential barrier, rotation cannot stop the inflow. The relativistic 
nature of a black hole gravitational field can be approximately described by 
the pseudo-Newtonian potential 92 = GM/(r — r g ) (Paczyhski & Wiita 1980) 
which is a good approximation for a Schwarzschild black hole. Similar poten- 
tials have been proposed to emulate the gravitational field of a rotating black 
hole (Artemova, Bjornsson & Novikov 1996). The pseudo-Newtonian models 
were studied by Paczyhski &; Bisnovatyi-Kogan (1981) and in many subsequent 
works. Recently, fully relativistic ID equations have been derived (Lasota 1994; 
Abramowicz et al. 1996). We here summarize the relativistic equations with 
some corrections (Beloborodov, Abramowicz &: Novikov 1997). 

First, we summarize notation. The Kerr geometry is described by the metric 
tensor gij in the Boyer-Lindquist coordinates x l = (t,r,6,<p) (given, e.g., in 
Misner et al. 1973). The four- velocity of the accreting gas is u l = (u l ,u r ,u 6 ,11?) 
with u 6 = assumed in the ID model. The angular velocity of the gas is 
£1 = uf /u* and its Lorentz factor is 7 = ^{—g 11 )^ 1 / 2 measured in the frame 
of local observers having zero angular momentum. The vertically integrated 
parameters of the disk are: surface rest mass density S, surface energy density 
U = Sc 2 +II (IT being the internal energy), and the vertically integrated pressure 
P. The dimensionless specific enthalpy is n = (U + P)/Sc 2 . F + is the rate of 
viscous heating, and F~ is the local radiation flux from the two faces of the disk. 
E, U, II, P, and F^ 1 are measured in the local comoving frame. 
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The main equations of a steady disk express the general conservation laws 
discussed by Novikov & Thorne (1973) and Page & Thorne (1974) in the context 
of the standard model. For advective disks, the assumption Q, = is replaced 
by the radial momentum equation, and the assumption F + = F~ is replaced by 
the energy equation including the advection term. 



Mass conservation: 



2TTrcu r T, 



-M 



Angular momentum: 



dr 



~2^~ 



2vY>ro r 



F~ 



ru v , 



where a r v = ^ g" ' g w \j '~9 a ^ 3 ^ is the shear 



Energy: 



p+ _ p- 



cu 



du n + p ds 



dr 



dr 



where F + = vY,/j,c 2 g rr g w (-5") 



Radial momentum: 



2 dr ; 2 dr 



5 V (O-Ok) 



1 dP F+u r 



c 2 T,fi dr c 3 T,fi ' 

where 0,^ are Keplerian angular velocities for the co-rotating (+) and counter- 
rotating (— ) orbits 



± 



r{2r/r g ) 1 / 2 ±a*r g /2' 



Here a* < 1 is the spin parameter of the black hole. 

The set of disk structure equations gets closed when the viscosity v, the 
equation of state P(H), and the radiative cooling F~ are specified. A standard 
prescription for viscosity is u = ac s H, where a is a constant and c s = c(P/U) l l 2 
is the isothermal sound speed. The half-thickness of the disk, H, should be 
estimated from the vertical balance condition. Near the black hole, the tidal 
force compressing the disk in vertical direction depends on 0, (e.g., Abramowicz, 
Lanza & Percival 1997). For Q ~ which is a good approximation for the 
vertical balance, one gets (e.g., Riffert & Herold 1995) 



H 2 = 



P 2r 3 

lj—y 



where 



J{a*,r) = 
5 



2(r 2 - a*r gy /2r~r~ + 0.75<r|) 
2r 2 — 3r 9 r + a*r gy /2r g r 



is a relativistic correction factor becoming unity at r > r 9 . 

The disk luminosity is related to M by (Beloborodov et al. 1997) 

OO / 

2-7T f I U m 

L~ = / u t F~rdr « Mc 2 1 + fi m — 

r s ^ 

where the index "in" refers to the inner transonic edge of the disk, in particular, 
fi m is the dimensionless specific enpthalpy at the inner edge. The radiative 
efficiency of the disk equals 

,^l-,,„(l-|). (8) 

Here, 6; n = c 2 + n™c is the specific binding energy at the inner edge. In the 
standard model, /j,{ n = 1 and rj = 6j n /c 2 . In the advective limit 77 — >■ 0, and hence 
A*in - ► (1 — frin/c 2 ) X - Note that one should not assume b in = for advective 
disks. In fact, the position and binding energy of the inner edge can be found 
only by integrating the disk structure equations. The difference fi{ n — 1 (which 
describes the relativistic increase of gas inertia due to stored heat) adjusts to 
keep rj rs 0. 



4. Super-critical disks 

Large accretion rates are expected in the brightest objects, such as quasars 
or transient GBHs during outbursts. When the accretion rate approaches the 
critical value m cr corresponding to the Eddington luminosity, the standard model 
becomes inconsistent. Inside some radius r t , the produced radiation is trapped 
by the flow and advected into the black hole, as the inflow time-scale here is 
less than the time-scale of photon escape from the disk, ~ tH/c, (Begelman & 
Meier 1982). The trapping radius can be estimated from the standard model of 
Shakura & Sunyaev (1973) by comparing the radial flux of internal energy 3c 2 M 
with the total flux of radiation emitted outside r. One then gets rt ~ 7r g for 
tin = rh cr = 12. The relative height of the standard disk, H/r, equals m/27 at 
the maximum. Advection thus starts to become important before the accretion 
flow becomes quasi-spherical, and it may be approximately treated retaining the 
vertically-integrated approximation. 

Advection reduces the vertical radiation flux F~ as compared to the stan- 
dard model and, as a result, the disk thickness stays moderate at super-critical 
accretion rates (Abramowicz et al. 1988; Chen & Taam 1993). This made pos- 
sible an extension of the ID model to the super-Eddington advection dominated 
regime, called "slim" accretion disk. The slim model may apply at moderately 
super-critical accretion rates. In the limit m 3> m CT there appears an extended 
region where the accretion flow has a positive Bernoulli constant, and the bulk 
of supplied gas may be pushed away by the radiation pressure to form a wind 
(Shakura & Sunyaev 1973). The behavior of the flow atm» rh cr still remains 
an open issue. Detailed 2D hydrodynamical simulations might clarify whether 
gas mainly falls into the black hole or flows out (Eggum, Coroniti & Katz 1988; 
Igumenshchev & Abramowicz 1998). 
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Figure 1. The temperature in the innermost region of the disk 
versus m for Schwarzschild (a* = 0, rh cr = 17.5, left panel) and Kerr 
(a* = 0.998, m cr = 3.11, right panel) black holes in the cases of a = 0.03 
(marked 1) and a = 0.3 (marked 2). A black hole mass M = 1O 8 M 
is assumed. The dotted curves show the blackbody temperature T e ff. 
(From Beloborodov 1998a.) 



The relativistic slim disk has been calculated in Beloborodov (1998a) for 
both non-rotating (a* = 0) and rapidly rotating (a* = 0.998) black holes. Like its 
pseudo-Newtonian counterpart, the relativistic slim disk has a moderate height 
up to m ~ 10rh cr due to advection of the trapped radiation. An important issue 
is the temperature of the accreting gas as the temperature determines the emis- 
sion spectrum. The simplest way to evaluate the temperature is by assuming 
that the gas is in thermodynamic equilibrium with the radiation density in the 
disk, w. Then T equals T e g = (w/a r ) l / A where, a r is the radiation constant. 
This was usually adopted in models of super-critical disks. The blackbody ap- 
proximation, however, fails when a > 0.03 (Beloborodov 1998a). The gas then 
accretes so fast and has so low density that it is unable to reprocess the released 
energy into Planckian radiation. As a result, gas overheats so that T S> T e g . 

The possibility of the overheating in near-critical disks was pointed out by 
Shakura & Sunyaev (1973) and discussed later (e.g., Liang & Wandel 1991; 
Bjornsson et al. 1996). In the standard model, the disk density scales as 
n oc m~ 2 , while the heating rate F + oc m. The resulting temperature of the 
overheated gas increases with m. The slim disk model allows one to follow this 
tendency toward the super-critical regime. Figure 1 shows the results for a mas- 
sive black hole in an AGN. The strongest overheating occurs at m ~ 3 — 4 m cr . 
At m > m cr the density increases, nam, and the temperature falls down. A 
similar overheating occurs in GBHs (M ~ 1OM ) if a > 0.1. 

The main process cooling the plasma in overheated disks with large M is the 
saturated Comptonization of bremsstrahlung photons (cf. Rybicki & Lightman 
1979). The plasma temperature is determined by the heating=cooling balance 



(Beloborodov 1998a), 



w = 



F+_ 
2H 



w s 



w 



w \ 



+ A\1- 

w p \ V w p i J 




(9) 



" 27 n 2 v / r is the free-free cooling rate, w p \ 



i r T 4 , and A is 
w p \ 10) this 



where tug = 1.6 • 10 

the Compton amplification factor. In the optically thick limit (w p \ ~ w 
equation yields w + = wg(l — w/w p \), while in the overheated case (w p \ S> w) it 
transforms into w + = Awg. The Compton amplification factor achieves A ~ 300 
in the hottest models. The efficient Compton cooling prevents a transition to 
an extremely hot two-temperature state (see also Bjornsson et al. 1996). 

The overheating will strongly affect the disk spectrum: the spectrum of a 
relativistic disk with rh ^ m cr can extend to the hard X-ray band, in contrast 
to the spectra of blackbody pseudo-Newtonian models (Szuszkiewicz, Malkan & 
Abramowicz 1996). 

Many GBHs and AGNs, however, accrete at sub-Eddington rates and their 
hard X-ray emission implies a different mechanism. 



5. Two-temperature disks 
5.1. The SLE model 

To explain the hard spectrum of the classical black hole candidate Cyg X-l, 
Shapiro, Lightman & Eardley (1976, hereafter SLE) suggested a two-temperature 
model of an accretion disk with protons much hotter than electrons, T p 3> T e . 
The hot proton component then dominates the pressure and keeps the disk geo- 
metrically thick. This in turn leads to a low density which scales oc H~ s for given 
rh and a. A low density implies a low rate of Coulomb energy exchange between 
the protons and the electrons and allows one to keep T p » T e self-consistently. 

The four main equations of the disk structure are the same as in the stan- 
dard model (see Section 2). The SLE disk has two temperatures. Therefore, the 
heating=cooling balance is now described by two equations: 

i) Energy balance for the protons. The model assumes that viscous dis- 
sipation heats preferentially the protons and then the bulk of this energy is 
transferred to the electrons through Coulomb collisions, i.e., 

F+ nkT p T e 3/2 

— — - = , where t ev « 17 

2H t ep ' p n 

is the Coulomb cooling time (with a Coulomb logarithm ~ 15, see Spitzer 1962). 
This equation yields 

Ve«10- 2 ^(£f (10) 

where 9 e = kT e /m e c 2 and 6 P = kT p /m p c 2 . 

ii) Energy balance for the electrons. The electrons cool mainly by upscat- 
tering soft radiation coming from an outer cold disk or from dense cloudlets 
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embedded in the hot flow. The cooling proceeds in the regime of unsaturated 
Comptonization (cf. SLE). In this regime, the disk parameters adjust so that 
y = 4# e max(T, r 2 ) ps 1. Typically, 6> e ;$ 1 and r ~ 1. The condition y « 1 gives 

^^fer-^ - i^r— • <"> 

Equations (10-11) yield an electron temperature that weakly depends on 
radius, 6 e ~ 0.1(am) _1 / 6 (r/r fl ) 1 / 4 . The complete set of equations allows one 
to express the other disk parameters as functions of radius (see SLE). In par- 
ticular, the proton temperature and the disk height are (H/r) 2 ~ 9 p /9 v i r ~ 
0.1a -4//3 m 2//3 (r/r 9 ) -1//4 where 6 v - lT ~ r g /r is the virial temperature. 

The model was further developed to include effects of pair production (e.g., 
Liang 1979b; Bjornsson & Svensson 1992; Kusunose & Mineshige 1992), a 
non-thermal particle distribution (Kusunose & Mineshige 1995), and cyclosyn- 
chrotron radiation (Kusunose & Zdziarski 1994). A relativistic version of the 
SLE disk around a Kerr black hole was calculated by Bjornsson (1995). 

The SLE model is in agreement with observed spectra of GBHs and AGNs. 
The shortcoming of the model is that it is thermally unstable (Pringle 1976; 
Piran 1978): the assumed energy balance would be destroyed by a small per- 
turbation of the proton temperature. An increase in T p would result in disk 
expansion in the vertical direction. Then the Coulomb cooling is reduced (and 
the heating rate is increased) and T p increases further, leading to an instability. 

5.2. The ADAF model 

In the SLE model only a fraction of the dissipated energy remains in the proton 
component to keep it hot, and the bulk of the energy is passed to the electrons 
and radiated. Ichimaru (1977) considered a disk of so low density that the 
heated protons are unable to pass their energy to the electrons on the time-scale 
of accretion. The protons then accumulate the energy and advect it. Such a 
flow is thermally stable. It emits only a fraction of the dissipated energy that 
has been passed directly to the electrons. 

If viscous dissipation heats mostly the protons, then the radiative efficiency 
of the accretion flow is very low, and the energy is either advected into the 
black hole (Ichimaru 1977; Rees et al. 1982) or transported outward in the 
form of a hot wind (Blandford & Begelman 1998, see also Begelman, these 
proceedings). The former case (ADAF) may be modeled in ID approximation 
(see Narayan et al. 1998 for a review). The latter case requires detailed 2D 
simulations (Igumenshchev & Abramowicz 1998) which are much more difficult: 
the very formulation of a 2D problem implies additional assumptions (local 
viscosity prescription, local heating rate, and 2D boundary conditions). 

The assumption that protons are heated preferentially has recently been 
assessed (see Quataert & Gruzinov 1998; Quataert, these proceedings, and ref- 
erences therein) . It was argued that the dissipation of Alfvenic turbulence in the 
disk may heat mostly protons only if the magnetic energy density is 10% of 
the proton pressure. If the magnetic field is the source of viscosity, one expects 
a < 0.1 in a two-temperature disk. 
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The necessary condition for ADAF models is that the time-scale for Coulomb 
cooling, t ep , exceeds the accretion time-scale, t a rs (af^x) -1 - It requires 

m < 10 2 a 2 ^ 2 . 

The typical electron temperature in the calculated ADAF models is ~ 10 9 K 
(e.g., Nakamura et al. 1997; Esin et al. 1998), and the corresponding maximum 
accretion rate is 

m m ax ~ 10a 2 . (12) 

From hydrodynamical point of view, the two-temperature ADAFs are sim- 
ilar to the super-critical disks. Again, most of the dissipated energy is stored 
inside the disk and it swells up to a quasi-spherical shape. The radial pressure 
gradients and deviations from Keplerian rotation are dynamically important, 
and the full set of differential equations should be solved to obtain a solution for 
the disk structure (Chen, Abramowicz & Lasota 1997; Narayan, Kato & Honma 
1997). At large distances, r 3> r g , the ADAF parameters have a power-law 
dependence on radius (Narayan & Yi 1994). An approximate description of the 
advective disk was given by Abramowicz et al. (1995) who assumed that the 
disk rotates with a Keplerian velocity in the pseudo-Newtonian potential. Re- 
cently, relativistic solutions have been calculated for ADAFs in Kerr geometry 
(Abramowicz et al. 1996; Peitz <fc Appl 1997; Igumenshchev, Abramowicz & 
Novikov 1998; Popham h Gammie 1998) and the expected spectra have been 
discussed (Jaroszyhski & Kurpiewski 1997). The ADAF model is applied mainly 
to low luminosity objects with suspected black holes, such as Sgr A*, nuclei of 
elliptical galaxies, and X-ray novae in quiescent state (see Narayan et al. 1998 
and references therein). 

The possibility of a two-temperature accretion flow in Cyg X-l and other 
GBHs in the hard state has again been addressed by Narayan (1996), Esin et al. 
(1998), and Zdziarski (1998). Their basic picture of accretion is the same as in 
the SLE model: at some radius r tr , the standard "cold" disk undergoes a tran- 
sition to a hot two-temperature flow which emits hard X-rays by Comptonizing 
soft radiation in the unsaturated regime. In contrast to SLE, radial advection 
of heat is taken into account, and the parameters (m and a) are chosen so that 
advection and radiative cooling are comparable. Then the flow is a good emitter 
and at the same time it may be stabilized by advection. It implies that the ac- 
cretion rate is just at the upper limit (eq. [12]), and a ~ 0.3 has to be assumed, 
as typically m ~ 1 in GBHs. The optical depth of the flow near the black hole 
is approximately 



Combined with the typical electron temperature T e ~ 10 9 K, this yields y ~ 1, 
just what is needed to be in the regime of unsaturated Comptonization that is 
believed to generate the X-ray spectrum. 

An advantage of a hot disk model is its consistency with the observed weak 
X-ray reflection in GBHs in the hard state (Ebisawa et al. 1996; Gierlihski et al. 
1997; Zdziarski et al. 1998; Zycki, Done & Smith 1998). Observations suggest 
that cold gas reflecting the X-rays covers a modest solid angle ~ 0.3 x 2n as 
viewed from the X-ray source. This is consistent with the outer cold + inner 
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hot flow geometry (for alternative models see Section 6.4). Provided the two 
conditions Tn rZi TT^niax and t ~ 1 are satisfied, and the inner radius of the cold 
disk, rtr, is properly chosen, the advective model may approximately describe 
the hard state of GBHs. A small increase (~ 10 %) in rh would then result 
in a transition to the soft state, as the hot flow collapses to a standard thin 
disk (Esin et al. 1998). A similar small decrease in m would lead to transition 
to a low-efficiency ADAF which may be associated with the quiescent state in 
transient sources. 

There is then a question, why the hard state is so widespread if it requires 
the fine-tuning of both m and a? In particular, why the hard state is stable in 
Cyg X-l while its luminosity varies by a factor of ~ 2? Variations in luminosity 
are presumably due to variations in m and, according to the ADAF model, the 
change in m must switch off the hard state and cause a transition to the soft 
or quiescent state. The latter has never been observed in Cyg X-l. Moreover, 
the observed slope of the X-ray spectrum in the hard state is kept constant 
with varying luminosity (Gierlihski et al. 1997), which indicates a persistent 
y-parameter of the emitting plasma. 

One possible explanation is that the disk has m > m max , being composed of 
two phases, hot and cold. The thermal instability of a hot disk at m > m max may 
continue to produce a cold phase until vn in the hot phase is reduced to ?rimax* 
Then a balance between the cold and hot phases might keep both conditions 
m and r ~ 1 for the hot phase (Zdziarski 1998). Hot accretion disks with 

cold clumps of gas have recently been discussed in different contexts (Kuncic, 
Celotti & Rees 1997; Celotti & Rees, these proceedings; Krolik 1998; Krolik, 
these proceedings). The cold clumps may also play a role as a source of seed 
soft photons for Comptonization (e.g., Zdziarski et al. 1998). 

The dynamics of a cold phase produced in the hot disk is, however, a very 
complicated and unresolved issue. In particular, it is unclear whether the cold 
phase may accrete independently with its own radial velocity, or if it rather is 
"frozen" in the hot phase. The cold clumps moving through the hot medium 
should be violently unstable, and a model for dynamical equilibrium is needed 
for a continuously disrupted/renewed cold phase. If the clump life-time exceeds 
the time-scale for momentum exchange with other clumps and/or with the hot 
medium, then the clumps should settle down to the equatorial plane and form 
a cold thin disk. One then arrives at the disk-corona model. 

6. The disk-corona model 
6.1. Magnetic flares 

That the observed X-rays can be produced in a hot corona of a relatively cold 
accretion disk extending all the way to the black hole was suggested a long time 
ago (e.g., Bisnovatyi-Kogan & Blinnikov 1977; Liang 1979a). Most likely, a low 
density corona is heated by reconnecting magnetic loops emerging from the disk 
(Galeev, Rosner & Vaiana 1979, hereafter GRV). This implies that the corona 
is coupled to the disk by the magnetic field (for alternative models, where the 
corona accretes fast above the disk, see, e.g., Esin et al. 1998; Witt, Czerny h 
Zycki 1997; Czerny et al., these proceedings). 
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The usually exploited model for the corona formation is that of GRV. Ac- 
cording to the model, a seed magnetic field is exponentially amplified in the 
disk due to a combination of the differential Keplerian rotation and the turbu- 
lent convective motions. The amplification time-scale at a radius r is given by 
*G ~ r/3v c where v c is a convective velocity. GRV showed that inside luminous 
disks the field is not able to dissipate at the rate of amplification. Then buoyant 
magnetic loops elevate to the corona where the Alfvenic velocity is high and the 
magnetic field may dissipate quickly. 

The coronal heating by the GRV mechanism is, however, not sufficient to 
explain the hard state of GBHs (Beloborodov 1999). The rate of magnetic energy 
production per unit area of the disk equals Fb = IBwbI^g where Wb = B 2 /8ir 
is the average magnetic energy density in the disk at a radius r. One can compare 
Fb with the total dissipation rate, F + = 3t rip c s , to get 

Fb H 



Here = B^Br/Ait and we took into account that B^j ' B r = c s /v c in the GRV 
model. Hence, the GRV mechanism is able to dissipate only a small fraction 
~ H/r <C 1 of the total energy released in the disk. 

Recent simulations of MHD turbulence indicate that the magneto-rotational 
instability efficiently generates magnetic energy in the disk (see Balbus & Hawley 
1998). The instability operates on a Keplerian time-scale, which is ~ H/r times 
shorter than tQ, and it produces magnetic energy at a rate Fb ~ F + . A low 
dissipation rate inside the disk would lead to buoyant transport of generated 
magnetic loops to the corona (GRV). The hard state of an accretion disk may be 
explained as being due to such a "corona-dominated" dissipation. By contrast, 
in the soft state, the bulk of the energy is released inside the optically thick disk 
and the coronal activity is suppressed. 

In the hard state, the corona is the place to which magnetic stress driving 
accretion is transported and released. Conservation of angular momentum reads 
(see eq. [4]) 

2Ht rip = — $7x5*. 

For a standard radiation-pressure-dominated disk it gives t rif ~ m p cVLK/&i:- 
If a large fraction, (, of F + is dissipated above the cold disk, then the disk 
height is reduced by a factor (1 — Q (Svensson & Zdziarski 1994). It implies 
that t r(fi is increased by a factor (1 — C) -1 - The magnetic corona is expected to 
be inhomogeneous and the main dissipation probably occurs in localized blobs 
where the magnetic energy density wb is much larger than the average t rip . The 
accumulated magnetic stress in such a blob may suddenly be released on a time- 
scale to ~ 10rfe/c (the "discharge" time-scale, see Haardt, Maraschi &: Ghisellini 
1994) where r& is the blob size. This produces a compact flare of luminosity 
L ~ r%w B /t - 

In this picture of the hard state, the binding energy of spiraling gas is 
liberated in intense flares atop the accretion disk. There is no detailed model for 
the magnetic flare phenomenon despite the fact that magnetic flares have been 
studied for many years in the context of solar activity. In particular, it is unclear 
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whether the flaring plasma should be thermal or non-thermal. Observations 
therefore play a crucial role in developing a model. Observed X-ray spectra 
of black hole sources suggest that the bulk of emission comes from a thermal 
plasma with a typical temperature kT ~ 50 — 200 keV and Thomson optical 
depth tt ~ 0.5 — 2 (see Zdziarski et al. 1997; Poutanen 1998). The narrow 
dispersion of the inferred T and tt indicates the presence of a standard emission 
mechanism. 

One possible scenario has been developed assuming a large compactness pa- 
rameter of the flares I = La^ /rbm e 10 — 10 3 . Then the flare gets dominated 
by pairs created in 7 — 7 interactions (e.g., Svensson 1986). The pairs pro- 
duced tend to keep tt ~ 0.5 — 2 and kT ~ 50 — 200 keV, in excellent agreement 
with observations. Yet, it is also possible that the flares are dominated by a 
normal proton plasma with tt ~ 1. Improved data should help to distinguish 
between the models. In particular, detection of an annihilation feature in the 
X-ray spectra would help. 

6.2. Compton cooling 

A flaring blob cools mainly by upscattering soft photons. We will hereafter 
assume that soft radiation comes from the underlying disk of a temperature 
T s . The additional source of seed soft photons due to the cyclo-synchrotron 
emission in the blob is discussed, e.g., by Di Matteo, Celotti & Fabian (1997) and 
Wardzihski & Zdziarski (these proceedings) . When a soft photon of initial energy 
e s passes through the flare, it acquires on average an energy Ae s , where A is the 
Compton amplification factor. A may also be expressed as A = (Ldiss + L s )/L s 
where Ldiss is the power dissipated in the flare and L s is the intercepted soft 
luminosity. 

The produced X-rays have a power-law spectrum whose slope T depends 
on the relativistic y-parameter of the blob, y = 4(6> e + 4^)tt(tt + 1). Leaving 
aside the effects of the (unknown) geometry of the blob, one may evaluate T 
by modeling radiative transfer in the simplest one-zone approximation in terms 
of an escape probability, as done, e.g., in the code of Coppi (1992). We have 
calculated the photon spectral index F using Coppi's code (see Figure 2). Within 
a few percent T follows a power law 

r « \y-^. (13) 

This empirical relation is simpler than the approximation of Pozdnyakov, Sobol 
& Sunyaev (1979), r » 1 + [2/(6> e + 3) - log r T ]/ log(12^ + 250 e ). 

One may also evaluate T as a function of the Compton amplification factor. 
Then the result depends on T s which is typically a few x 10 6 K in GBHs and a 
few x 10 4 K in AGNs. The corresponding dependences T(A) are shown in Figure 
2b. To high accuracy (~ 3 — 4 %), the results can be approximated as 

T^ 7 -(A-l)- s , (14) 

where S « 1/6 for GBHs and S « 1/10 for AGNs. Formula (14) is more accurate 
than the estimate of Pietrini k Krolik (1995), r « 1 + l.Q{A - 1)~ 1/4 , where the 
dependence on T s is neglected. 
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0.5 1 0.5 1 1.5 2 

log y log (A-l) 

Figure 2. (a) X-ray spectral index, T, versus y-parameter. Open 
circles - kT s = 5 eV, kT = 50 keV, open squares - kT s = 5 eV, 
kT = 100 keV, filled circles - kT s = 200 eV, kT = 50 keV, filled squares 
- kT s = 200 eV, kT = 100 keV. The line shows the approximation (13). 
(b) T versus A — 1 = L^ ss /L s . The two lines show the approximation 
(14) with 5 = 1/10 and 5=1/6. 

6.3. The feedback of X-ray reprocessing by the disk 

The flares illuminate the underlying accretion disk that must reflect/reprocess 
the incident X-rays. The disk produces important features in the observed X-ray 
spectrum, such as the Fe Ka line and the Compton reflection bump, as exten- 
sively discussed in these proceedings. Even more important, the disk reprocesses 
a significant part of the X-ray luminosity into soft radiation. As the flare is ex- 
pected to dominate the local disk emission, the reprocessed radiation becomes 
the main source of soft photons. Then the flares are "self-regulating": the flare 
temperature adjusts to keep A = (L s / L)~ l where L s /L is a fraction of the flare 
luminosity that comes back as reprocessed radiation (Haardt & Maraschi 1993; 
Haardt et al. 1994). The resulting spectral slope is determined by the feedback 
of reprocessing, L s /L. Detailed calculations of the predicted X-ray spectrum 
were performed by Stern et al. (1995) and Poutanen & Svensson (1996) and 
applied to Seyfert 1 AGNs (see Svensson 1996 for a review). 

The calculated models are, however, in conflict with observations of Cyg 
X-l and similar black hole sources in the hard state (e.g., Gierlihski et al. 1997): 
i) The observed hard spectrum corresponds to a Compton amplification factor 
A £ 10 and implies soft photon starvation of the hot plasma, L s <C L. The 
model predicts A ^ 5 unless the active blobs are elevated above the disk at 
heights larger than the blob size (Svensson 1996). ii) The model with elevated 
blobs would yield a strong reflection component, R = VL/2it ps 1, where Vt is the 
solid angle covered by the cold matter as viewed from the X-ray source. The 
reported amount of reflection is small, R ~ 0.3. 
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The weak reflection and soft photon starvation may be explained if the cold 
reflector is disrupted near the black hole. This would agree with the idea that 
accretion proceeds as a hot two-temperature flow in the inner region as discussed 
in Section 5. One may fit observed spectra with a toy model of a hot central 
cloud upscattering soft photons supplied by a surrounding cold disk or by dense 
clumps inside the hot region (e.g., Poutanen, Krolik & Ryde 1997; Zdziarski et 
al. 1998). The transient soft state is then explained as being due to shrinking 
of the hot region, so that the cold disk extends all the way to the black hole and 
the bulk of the energy is dissipated inside the optically thick material of the disk 
(Poutanen et al. 1997; Esin et al. 1998). 

However, the weak reflection does not necessarily imply that the inner cold 
disk is disrupted in the hard state. One suggested alternative is that the ap- 
parent weakness of the reflection features is due to a high ionization of the 
upper layers of the disk (Ross, Fabian & Young 1998). Another alternative is 
that the emitting hot plasma has a bulk velocity directed away from the disk 
(Beloborodov 1999). Mildly relativistic bulk motion causes aberration reducing 
X-ray emission towards the disk. It in turn reduces the feedback of reprocessing 
and leads to a hard X-ray spectrum. 

The coupling between the flare and the underlying disk can be approxi- 
mately described assuming that the luminosity emitted downwards within an 
angle cos -1 fi s comes back to the flare. The effective fi s depends on the flare 
geometry. E.g., a slab geometry of the active region corresponds to ^ s = 0, and 
an active hemisphere atop the disk has /i s ~ 0.5. The feedback factor, L s /L, of 
a flare of luminosity L atop the disk is determined by three parameters: 

• The geometrical parameter fj, s . 

• The bulk velocity in the flare, (3 = v/c (assumed to be perpendicular to 
the disk). 

• The disk albedo, a. x = 1 — a represents the efficiency of reprocessing of 
the incident X-rays. 

In the static case (j3 = 0), the reflection R = 1 and L s = L%(1 — /i s )/2. 
The corresponding amplification factor is A = 2/^(1 — /U s ). With increasing 
/3, A increases, and R is reduced. The impact of a bulk velocity on A and R 
is summarized in Figure 3a for several fi s and system inclinations 9. In the 
calculations, we assumed a typical albedo a = 0.15 (e.g., Magdziarz & Zdziarski 
1995). Note that the observed R may be further reduced because the reflected 
radiation is partly upscattered by the blob. 

One can evaluate the spectral index of a flare using equation (14). For a 
typical fi s ~ 0.5 one gets T « 1.9B- 5 for GBHs and T « 2B' 3 for AGNs, 
where B = 7(1 + 0) is the aberration factor due to bulk motion (Beloborodov 
1999). E.g., for Cyg X-l in the hard state, both the spectral slope V ~ 1.6 and 
the amount of reflection R ~ 0.3 can be explained assuming (5 ~ 0.3. 

6.4. The ejection model 

The inferred (3 > implies that the flares are accompanied by plasma ejection 
from the active regions, in contrast to the static corona model. In fact, one 
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Figure 3. (a) Top: Compton amplification, A, of the reprocessed 
radiation as a function of (3. A disk albedo a = 0.15 is assumed. 
Bottom: Apparent amount of reflection, R. The dotted line shows 
expected A and R for j3 = 0.3. (From Beloborodov 1999.) 
(b) Best fits to T and R for a number of AGNs (open circles), Cyg 
X-l (filled squares), and GX 339-4 (filled circles). The data are from 
Zdziarski et al. (1999). The solid curve shows the correlation R(F) 
predicted by the ejection model for AGNs (assuming fi s = 0.55 and 
inclination 8 = 45°). R(T) predicted for AGNs by the hot disk model 
(see Zdziarski et al. 1999) is shown by the dotted line. 



should expect bulk motion of a flaring plasma on theoretical side, especially 
if the plasma is composed of light e ± pairs. There is at least one reason for 
bulk acceleration: the flare luminosity, L, is partly reflected from the disk, and 
hence the flaring plasma is immersed in an anisotropic radiation field. The net 
radiation flux, ~ L/r%, is directed away from the disk and it must accelerate 
the plasma. The transferred momentum per particle per light-crossing time, 
r^/c, is ~ lm e c where I = Lai/ri,m e c 3 is the compactness parameter of the 
flare. Hence, the acceleration time-scale is t a ~ I ^Tfo/c for a pair plasma and 
t a ~ (m p /m e )l~ 1 ri } /c for a normal proton plasma. The shortness of t a for a pair 
plasma implies that the pair bulk velocity saturates at some equilibrium value 
limited by the radiation drag. Using a simple toy model, one may estimate the 
expected velocity to be in the range f3 ~ 0.1 — 0.7 (Beloborodov 1999). A proton 
plasma may also be accelerated to relativistic velocities if the flare duration 
exceeds t a , which is quite probable. 

The very magnetic dissipation may be accompanied by pumping a net mo- 
mentum into the flare at a rate ~ L/c. When the stored magnetic energy gets 
released, the heated plasma may be ejected both toward and away from the disk. 
Again, the large compactness parameter implies efficient momentum transfer to 
the plasma, ~ lm e c per particle per light crossing time. Possible ejection toward 
the disk corresponds to < 0. 
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Plasma ejection from magnetic flares is likely to occur in both GBHs and 
AGNs. The plasma velocity may vary. An increase in (3 leads to decreasing 
R and T. A correlation between R and V is observed in GBHs and AGNs 
(Zdziarski, Lubihski & Smith 1999; Zdziarski, these proceedings) and it is well 
reproduced by the ejection model, see Figure 3b. The theoretical curve is plotted 
for AGNs, assuming a disk albedo a = 0.15, fi s = 0.55, and inclination 9 = 45°. 
These parameters are probably the most representative. The dispersion of data 
around the curve might be due to a dispersion in fi s and inclination angles. The 
velocity, (3, varies from — 0.2to0.75;i? = l corresponds to (3 = 0. 

The two model curves shown in Figure 3b are calculated for AGNs (<5 = 1/10 
in eq. [14]). For GBHs, V should be systematically smaller as a result of a higher 
energy of seed soft photons. This tendency is seen in Figure 3b: the Cyg X-l 
and GX 339-4 data are shifted to the left as compared to the AGN data. 

Outflows are usually expected in radio-loud objects where they have been 
considered as a possible reason for weak Compton reflection (Wozniak et al. 
1998). Figure 3b may indicate that plasma acceleration in X-ray coronas of 
accretion disks is a common phenomenon in black hole sources. Note also that 
fast outflows may manifest themselves in optical polarimetric observations of 
AGNs (Beloborodov 1998b; Beloborodov & Poutanen, these proceedings). 

7. Concluding remarks 

Disk-like accretion may proceed in various regimes. Any specific model is based 
on assumptions including prescriptions for the effective viscosity, the vertical 
distribution of the ion/electron heating, the energy transport out of the disk, 
etc. Large uncertainties in these prescriptions allow one to produce a lot of 
models, and observations should help to choose between them. Hard X-ray 
observations play an important role in this respect. They indicate that a large 
fraction of the accretion energy is released in a rarefied plasma (corona) where 
the X-rays are generated by unsaturated Comptonization of soft photons. In 
most cases, the accreting mass is expected to be concentrated in a "cold" phase, 
probably forming a thin disk embedded in the corona. Its gravitational energy 
transforms into the magnetic energy that is subsequently released in the corona. 

Most likely, the dissipation mechanism is related to highly non-linear MHD 
which usually produces inhomogeneous and variable dynamical systems. In par- 
ticular, the energy release in the corona is likely to proceed in magnetic flares 
generating observed temporal and spectral X-ray variability. Idealized hydro- 
dynamical models are unable to describe this process. The difficulty of the 
problem is illustrated by the activity of the Sun, where theoretical progress is 
modest despite the fact that detailed observations are available. In accretion 
disks, magnetic fields are likely to play a crucial dynamical role and various 
plasma instabilities may take place. The instabilities probably govern the dis- 
tribution of the plasma density and the heating rate. 

Given the difficulty of the accretion physics, studies of phenomena having 
specific, observationally testable implications, are especially useful. In partic- 
ular, the reflection and reprocessing of the corona emission by the cold disk 
provides diagnostics for accretion models. We here discussed a plausible phe- 
nomenon - bulk acceleration of the flaring plasma in the corona - that strongly 
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affects the reflection pattern and the radiative coupling between the corona and 
the cold disk. 
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